The n — > oo limit of 0(n) models on graphs 
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Thirty years ago, Stanley showed that an 0(n) spin model on a lattice tends to a spherical 
model as n — » oo. This means that at any temperature the corresponding free energies coincide. 
This fundamental result, providing the basis for more detailed studies of continuous symmetry spin 
models, is no longer valid on more general discrete structures lacking of translation invariance, i.e. 
on graphs. However only the singular parts of the free energies determine the critical behavior of the 
two statistical models. Here we show that such singular parts still coincide even on general graphs 
in the thermodynamic limit. This implies that the critical exponents of 0(n) models on graphs for 
n — ► oo tend to the spherical ones and therefore they only depend on the graph spectral dimension. 
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The study of spin models on non crystalline structures is an intriguing and complex problem in statistical me- 
. chanics. This is fundamentally due to the lack of translational invariance and of a natural definition for the system 
dimensionality The former gives rise mainly to technical difficulties, arising from the impossibility of using such 
t-H , a powerful tool as Fourier transforms. The latter involves some deeper questions concerning the role of large scale 
geometry in phase transitions on important real structures such as amorphous solids, glasses, polymers and fractals. 

Recently, these problems have been successfully addressed using graph theory. A graph, i.e. a network composed 
of sites and links connecting nearest-neighboring sites, is the most suitable geometrical model to describe an irregular 
| system consisting of spins coupled by local interactions. In particular, algebric graph theory gives very interesting 
results when dealing with continuous symmetry models, due to the evidence for deep relations between their critical 
QO ' behavior and the small eigenvalues spectrum of the Laplacian operator. 

OS . One of the most representative models in this class is the 0{n) model, which describes a classical n— dimensional 
spin vector of fixed length. It is not in general exactly solved, but a few analytic results have been obtained on 
lattices. For example the lower critical dimension for spontaneous symmetry breaking is known from the Mermin- 
Wagner theorem and the n — > oo limit, which corresponds to the spherical model, can be exactly solved on any lattice. 
These results provide the basis for a qualitative understanding of the phase diagram as well as for the 1 /n expansion 
of thermodynamical quantities and, in particular, of critical exponents. 

As for graphs, while an extension of the Mermin- Wagner theorem to non translation invariant structures has been 
proven 0], the n — > oo result for 0(n) models heavily relies on the invariance properties of the lattice and therefore 
cannot be easily generalized. Indeed the equivalence between an infinite components model with local constraints 
• i— i , and a model with a global constraint on the spin is lost as a consequence of the lack of symmetry. The extension of 
this result to a generic discrete structure would be an important step in the comprehension of continuous symmetry 
models and in particular of their phase diagrams and universality classes. In this letter we deal with such a problem. 

We show that on a graph, although the free energy of 0(n) models in the n — ► oo limit at a generic temperature 
is not in general equivalent to the corresponding free energy of the spherical model, the singular parts of the two 
free energies concide. Therefore on a generic discrete structure the equivalence is retrieved asymptotically near the 
critical point and it holds for critical exponents, which can be shown to depend only on the spectral dimension of the 
network. This result strongly supports the possibility of a global geometrical characterization of irregular networks 
in the critical region, affecting all phenomena related to large scales. 

The ferromagnetic classical 0(n) Heisenberg model is defined by the Boltzmann weight exp(— f3H n ), where 

//,/S ^./,.S,-S, h-^S,. (I) 
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the sum extends to all links of a certain graph Q with N sites, Jy > are ferromagnetic interactions, which may vary 
from link to link, and Si is an n— dimensional vector of fixed length normalized by Si • Si = n. The Hamiltonian (|l|) 
is equivalent to 

H n [S] = i £ ~ S ^) 2 ~ h ' E S *< ( 2 ) 

i,j i 

apart from an additive constant. In the following we will use this form for H n which simplifies the mathematical 
proof. The free energy per component is defined by 

where the partition function reads 

z = I n s t s * - s *- n ) ds * e ~ m - w 

J i 

We shall assume the existence of the thermodynamic limit N — > oo which turns Q into an infinite graph within a 
certain class to be better specified in the sequel. 

The classical fundamental result originally due to j| applies when Q is a regular lattice (e.g. a hypercubic lattice); 
it establishes a rigorous relation between 0(n) models in the n — > oo limit and the spherical model. The latter is 
defined, on the same graph Q of the 0(n) model, by the Boltzmann weight e~^ H with Gaussian Hamiltonian 

hS = \ E Mb -to? (5) 

<ij> i 

and the spherical constraint J^. <f> 2 — N, where <f>i are real scalar variables. In the thermodynamic limit N — > oo we 
are allowed to replace (||) by 



H s 
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\ E J v(& - + \ m2 - 1) - & (6) 



where the "mass" parameter m 2 is fixed to be a precise function of temperature and magnetic field, i.e. m? = fx(/3, h), 
by the spherical constraint "on the average" 
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Since H s is Gaussian, one finds (here and in the following we restrict to the case h = 0) 



f(™ 2 ) =~^ 00 ^p l0 % zS = ^logdet(L + m 2 ) - ?L , (<^.) = l( L + m a)-i ( 8 ) 

where -Ly- = Ji<5y — Jy, with Ji = J^. Jy, is a discrete Laplacian operator on CJ. At this stage it is convenient to better 
specify the class of infinite graphs and couplings with which we are concerned: we shall assume that Q can be naturally 
embedded in a finite dimensional Euclidean space and that Ji is uniformely bounded over Q. For our class of graphs a 
generic quantity qt, related to a single point i, can be averaged in a unique way over Q by [q]g = Unizv— >oo l/iV^. 
The measure \Q'\ of a subgraph Q' is given by the average value of its characteristic function \f > defined by xf = 1 
if i S Q' and Xi = otherwise: = [x]e- The asymptotic behavior of the model in the massless limit m 2 — > is 
related to the spectrum of the Laplacian operator at low eigenvalues Q . In particular we have the following singular 
behaviors as m 2 — > 



1 N 

sing [(£ + m 2 )-^ = sing lim - £(L + m 2 )^ = C" (m 2 )^ 1 (9) 



( = 1 



which define the spectral dimension d of the graph. Together with eqs. ([?]) (which nows reads [(4> 2 )]g = 1) and (||), 
eq. (|h determines the the asymptotic form of the function m 2 — n((3, 0) = £i(/3), that is /j,(/3) ~ /3 2 /^ _2 \ /3 — > oo for 
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d < 2 and ft(P) ~ (fl c - j3) 2 / {3 - 2 \ $ ->• #T for 2 < d < 4 and /i(/3) ~ & - /3 for d > 4, where /3 C = [L _1 ] g . In turns 
this implies that the free energy of the spherical model (see eq. (|8|)) has the asymptotic form near the critical point 

/ s ( M (/3))^/ s (0) + O(M/3))- (10) 

It is important to observe that for most infinite graphs, if Q' is a subgraph of Q with non-zero measure, \Q'\ > 0, 
the infrared singularities of [(L + m 2 ) 1 ]^ and [(L + m 2 ) 1 ]gi concide and are described by the same spectral dimension 
d and the same coefficient C of the singular part || . We call these graphs "pure graphs" . This is the case for all 
fractals (such as Sierpinski gasket and carpet, t-fractals and so on), bundled graphs (comb lattices, brushes, ...) and 
many others (e.g. NTjj 0). However some peculiar cases of macroscopically inhomogeneous networks, obtained for 
example by sticking togheter two pure graphs with different spectral dimensions by a zero-measure set of links, can 
present different singular behaviors on different non zero measure subgraphs. Now, it can be shown J3] that in these 
cases, which we will call "mixed", the number of such subgraphs must be finite and the Gaussian free energy for 
the whole graph is simply the weighted sum of the Gaussian free energies of the subgraphs. This reduces the global 
problem to the study of a finite number of pure problems. Therefore, in the following we will restrict to the pure case. 

Let us now come back to the equivalence between the O(oo) model and the spherical one. The proof of the classical 
result mentioned above deeply relies on the translation invariance of lattices with simple elementary cells whenever the 
couplings Jij are constant over the lattice. It holds in the thermodynamic limit and uses a saddle points technique j^j 
that cannot be generalized to a generic discrete structure. Here, we will follow a more general and flexible approach to 
the infinite component limit. In fact it can be proven that, provided < e < Jij < Jm < °o (bounded ferromagnetic 
couplings) §]: 

- i- log[(m 2 + J M )] < /„ - f s (m 2 ) < K^YKtf + V 

i 

where the constants K\ and do not depend on m 2 , N and n, f n is the free energy per component for the 0(n) 
model defined by M), f s (m 2 ) is that defined in eq. (||) and the average in the right has Boltzmann weight e _,9ffS , with 
the Hamiltonian (g) with generic m 2 , i.e where m 2 is not yet fixed by the average spherical constraint (Q). Notice that 
inequality (|ll| ) contains as a particular case the hypercubic lattice result. Indeed on a translation invariant structure 
((f> 2 ) does not vary from site to site and the spherical constraint (|7|), corresponding to m 2 = /i(/3), allows the right 
hand side of ( jll| ) to vanish for n — > oo, implying the coincidence of the two free energies. 

In the case of a generic discrete structure (</> 2 ) does change from site to site and one would in general need site- 
dependent squared masses mf to enforce (0 2 ) = 1 for every i and obtain an analogous coincidence of the two free 
energies in the infinite components limit. Thus in the thermodynamic limit one would have to solve an infinite 
number of equations which could not be reduced to a single one due to the lack of translation invariance, forbidding 
the equivalence with the spherical model where only the average constraint (Q) holds and is solved through a single 
global parametrization m — /i(/3). 

The key point of our approach is that we are interested in the critical properties of the two models. Therefore we 
will require from the beginning only the concidence of the singular parts of the two free energies. This requirement has 
two main consequences: the constraint equations are replaced by "constraint inequalities" and, above all, the original 
models can be replaced by modified models with the same singular parts of the free energies, therefore belonging to the 
same universality classes. Then we will show that a particular global choice m 2 = p,(/3) corresponding to the solution 
of the constraint equation for a rescaled spherical model on Q do satisfy the constraint inequalities and therefore gives 
the critical behavior of corresponding rescaled 0(n) model in the n — > oo limit. Finally, the rescaled models will be 
shown to belong to the same universality class of the original ones. 

The first step of our proof consists in obtaining the constraint inequalities. 

Let us consider the n — > oo limit of eq. (Jill) , which we write as 



and suppose we could show that 



f S (m 2 ) < /oo < f s (m 2 ) + K^ 2 - 1)% (12) 



[<0 2 - l)\ < o(p((3)) (13) 

near the critical point of the spherical model. Then comparing with the asymptotic form ( |To| ) of f s (n((3) itself would 
immediatly prove that the two free energies have both the same value at the critical point and the same singular parts 
near it. 

Now, for the class of infinte graphs under consideration it can be shown that (|l3|) follows from the infinite set of 
linearized inequalities 
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m 2 -i)»h>\<o{m) (14) 

where the average is taken over every subgraph Q> with \Q'\ > 0. We shall call (0) the "constraint inequalities" by 
analogy with the usual approach. 



As a matter of fact, it is in general impossible that inequalities (14) are verified by the solution /z(/3) of the global 



constraint for the standard spherical model on Q. However, we prove the following: 

• it is possible to find a modified set of couplings {JL} and site-dependent masses {mf}, both functions of m 2 , 
which define two modified models on Q with Hamiltonians H and H' n such that the corresponding constraint 
inequalities (|l4|) are satisfied by the solution of the (modified) spherical model on G, namely m 2 = ft(ft), with 
p.(P) = 0(/i(/3) near criticality. Therefore, from the inequality (|lj), the two corresponding free energies f^iP) 
and f' S (fl(P)) have the same singular part; 

• the rescalings {Jij} — ► {J'ij} an d {™ 2 } — ► {mf} do not affect the singular parts of the free energies for both 
the 0(n) and the spherical models and therefore such modified models belong to the same universality classes 
as the original ones. 

To prove the first point let us introduce the set of couplings and masses 

Jij = Jij\J[ai + bi(m 2 )][a,j +bj(m 2 )] , mf = m 2 [ai + bi(m 2 )] (15) 

where 0<e<ai<A<cx) and the bi(m 2 ) are functions to be determined which vanish for m 2 — > 0. Such a 
rescaling can be seen as a bounded rescaling of the scalar spin variables fa — > 0^ = fa[cii + bi(m 2 )] with the original 
Hamiltonian (0). Therefore, the inequalities ( |l4"|) for the modified model H' s can be written in terms of the (/)[. 
Now setting mr = the constants aj and the functions bi(fi) can always be chosen self-consistently to satisfy the 
inequalities to o(/i(/3). This can be done at the same time for every subgraph Q' with \Q'\ > 0, as required by eq. 
(|l4|), because the leading singularity of [(0 2 )]g', including its coefficient, is universal, that is equal to that of [((t> 2 )]g- 
Then for n — » oo the 0(n) model with coupling and the spherical model, defined by the same set of J[j satisfy 
the constraint inequalities. Therefore, their free energies has the same singular part and the two models have the 
same critical properties. 

Now, let us consider the universality classes of these models. The critical behavior of the spherical model on a 
graph can be expressed in term of its spectral dimension d, defined by (|^) . This has been shown to be invariant under 
bounded rescaling of couplings Jy such as the ones corresponding to the set J^. Therefore, the two spherical 
models described by the and JL belong to the same universality class. On the other hand by the generalized 
Griffith inequalities Q the correlation functions of the 0(n) model corresponding H' n has the same critical behavior 
of the correlation function of the model H n . Then they belong to the same universality classes. This completes our 
proof. 

Our result on the equivalence of critical regimes for the spherical and O(oo) models on graphs provides a key 
step for the comprehension of phase transitions on general networks. The spherical critical exponents, which are 
exactly known and depend on the spectral dimension of the graph, are the starting point for the determination of 
the corresponding ones for 0{n) models, via the 1/n expansion. This expansion is a fundamental tool to analyze the 
existence of geometrical universality classes on graphs: indeed if its coefficients can be shown to depend only on the 
spectral dimension d, then their existence is proved. 
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